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I . INTRODUCTION 



A . GENERAL 

Numerical path planning has been studied quite a bit in 
the last few years. One problem of numerical path planning 
involves finding the optimal path from a given starting point 
to a goal point through a plane that has been subdivided into 
weighted regions . This problem is known as the Weighted Region 
Least Cost Path (WRLCP) problem. The best path can be 
minimizing some cost (e.g. energy or time) . 

Several approaches can be used to tackle this WRLCP 
problem. Some of these approaches and their pros and cons are 
presented in chapter II. 

One of the most significant constraints has been the 
ability to evaluate in real time the optimal path through a 
weighted region. Traditional sequential algorithms can 
quickly become overloaded with the sheer number of 
calculations required for a realistic problem. Consequently, 
several algorithms for solving this type of problem have been 
researched in recent years. These algorithms have, to a 
varying degree, potential for parallelization. Parallel 
algorithms have a tremendous advantage in speed, but also may 
open the door to new problems, especially in mapping, 
scheduling and coordinating the different parallel 
activities. Multicomputer networks are no exception. The 
principal goal in parallelizing an algorithm is to sustain a 
close-to-linear speedup in processing time. This is no trivial 
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task, since the processor communication reduces the speedup. 
This, and the fact that parallelizing an algorithm may 
increase its overall processing time (summed over the total 
network of processors) make this problem challenging. 

We investigated in this thesis some of the basic approaches 
with their potential for parallelization, and implemented one 
of these algorithms, first sequentially and then in parallel 
on a Transputer network. 



B. THESIS OUTLINE 

Chapter II introduces an analysis of the different 
parallel path approaches which have been investigated 
recently. The algorithm which appears most promising for 
solving the weighted region least cost path problem, which is 
a grid algorithm, has been selected for further study. In 
chapter III, implementation and some modifications of this 
algorithm is presented with an illustrative example. Chapter 
IV provides a description of the Transputer network and ADA. 
In chapter V, this approach has been developed into a parallel 
algorithm and implemented on an INMOS Transputer network using 
ADA. The algorithms ability to use parallel computing to solve 
arbitrary WRLCP problems has been investigated. Finally, 
these results will be compared with the results of single 
processor, and with previous results on a multicomputer 
network. Conclusions and recommendations for further research 
are offered in chapter VI. 
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II . BACKGROUND 



We overview in this chapter some of the basic approaches 
for solving the WRLCP problem and we investigate their 
potential for parallelization. 

A. SINGLE SOURCE SHORTEST PATHS (SSSP) 

1. Dijkstra's Algorithm 

The algorithm for SSSP presented in Figure 1, which was 
developed by Dijkstra in 1959, is the starting point for this 
problem [MAN89] . 



Algorithm Single_Source_Shortest_Pathts (G,v); 

Input : G=(V.E) (a weighted directed graph) , and v (the source vertex). 
Output : for each vertex w, w.SP is the length of the shortest path from v to w. 
(afl lengths are assumed to be nonnegative. ) 

begin 

for all vertices w do 
w.mark := false; 
w.SP := infinite; 
end loop; 

v. SP := 0; 

while there exists an unmarked vertex do 

let w be an unmarked vertex such that w.SP is minimal; 

w. mark := true; 

for all edges (w,z) such that z is unmarked do 
if w.SP + lengdt(w,z) < z.SP then 
z.SP := w.SP + length(w,z) 
end if; 
end loop; 
en ejid loop; 



Figure 1: Dijkstra's Algorithm 



In Figure 2, a small example is presented to 
demonstrate the algorithm. The first line includes only paths 
of one edge from v (the source) . The shortest path is chosen, 
in this case, leading to vertex a. The second line shows the 
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update of the paths including now all paths of length one from 
either v or a, and the shortest path now leads to c. A new 
vertex is added in each line, and the current known shortest 
paths from v are listed to every vertex. The underlined 
distances are those that are known to be the shortest . The 
algorithm keeps adding new vertices to the selected list until 
all vertices are added. 
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Figure 2: An Example Of The SSSP Algorithm 
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The algorithm can be easily extended from directed 
to undirected form. 

2. Parallelization Of SSSP 

Dijkstra's algorithm has been parallelized by Moore 
[QUI87], He devised two parallel algorithms. The first 
algorithm makes the for loop (in Figure 1) parallel, which 
explores the outgoing edges from a given vertex, the second 
method is to parallelize the while loop (in Figure 1); that 
is, at any one time in the execution of the algorithm there 
are probably many vertices in the queue. The parallelizability 
of the first method is restricted by the number of edges 
outgoing from each vertex. On the other hand, the second 
method performs larger tasks, that produces good speedup. More 
detailed information can be found in the given reference. 

In the parallel algorithm based on the second method, 
a queue is used. That queue is initialized with the source 
point, and then a number of asynchronous processes are 
created. Each of these processes goes through the steps of 
deleting a node from the queue, examining its outgoing edges, 
and inserting into the queue the nodes to which shorter paths 
have been found. In Figure 3 an example is presented in which 
the number of nodes in the queue shows actually the number of 
processes that can be parallelized. Distances are kept until 
they are reached from another direction. When there are more 
than one edge reaching to a node, minimum cost is chosen. 
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Figure 3: Parallelized SSSP Algorithm Example 



B. RAY TRACING 

1. Snell's Law of Refraction 

Snell's law [HEC87] defines the path a ray of light 
takes as it passes from one medium to another. The ray is 
refracted across the border between the media according to the 
following equation, 

rtjSinOj = n 2 sinQ^ 

where n 2 are the indices of refraction and the angles are 
of incidence and refraction respectively (Figure 4) . Fermat's 
principle which implies Snell's law states that light follows 
a path between two points such that it takes the minimum time. 
It can be proven (for example see [RIC87]), that a similar 
principle govern the path a particle takes across two regions, 
in which the speed of the particle is uniform in each region. 
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Therefore, we can apply the principles of optics to solve the 
WRLCP problem by assuming regions as optical media and weights 
as indices of refraction. 




Figure 4: Snell's Law Of Refraction 

2. Implementation 

The main step of this implementation is to shoot a ray 
from the source and see where it is going to land. We keep 
doing this for initial rays with different angles until some 
of the rays hit the goal. In essence, the boundaries border 
the weighted regions, the index of refraction depends on the 
weights and the rays refract on the border to border until it 
reaches the goal. At every boundary of two regions the ray 
obeys Snell's law of refraction. Among those rays that hit the 
goal point we obtain the WRLCP. In Figure 5, at the source 
point the short rays are just to give an idea about different 
set of angles and continuing three rays are given to make a 
clear example. 
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Figure 5: Ray Tracing 

Through different weighted regions, as the path Sklmnop 
misses the goal, the minimal paths SabcdeG and SfghijG hit it, 
and one of them is chosen as WRLCP after comparison. 

3. Parallelization 

The Snell's law based algorithm can be easily 
parallelized by assigning a different set of angles to each 
processor. If a solution is found, it is optimal. The cost of 
finding the intersection point of a ray with a triangle is not 
large. However, the algorithm suffers several drawbacks: The 
ray inversion, possible presence of blind regions, and use of 
expensive trigonometric functions. More details can be found 
in [RIC87] . 
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C. EDGE SLICING 

1. Implementation 

In this technique, every edge is divided into a number 
of segments of equal length. Some points are placed 
equidistantly on every map edge in the triangulated plane. A 
graph is constructed by connecting every two points on two 
edges belonging in the same triangle. The distance between two 
consecutive points is made proportional to some function of 
the costs of the two regions separated by that edge. A graph 
is constructed whose nodes are the original triangles' 
vertices plus the points which divides the edges into 
segments. The edges of this graph are the original triangles' 
edges plus the lines connecting every two non-colinear points 
in the same graph. In Figure 6, edge slicing is shown for the 
case where an edge is divided into two equal segments. 

2. Parallelization 

After the edges in a graph are sliced, the WRLCP problem 
is reduced to a SSSP problem, and the technique used for 
parallelizing SSSP can be invoked. A variant of this algorithm 
is shown in [KIN91] . 

D. GRID ALGORITHMS 

Another approach for approximating the WRLCP problem is to 
model the terrain as a grid. Simply, a grid is laid over the 
terrain. The map is divided into equidistant grid points. The 
weights in the regions are transferred as edge costs. Figure 
7 illustrates how a node is connected to different number of 
neighbors (we implemented a node with 4-neighbors, but it is 
easy to extend the implementation for more neighbors) . 



9 




Figure 6: Edge Slicing (2-Point) 



Finally, the shortest path analysis is performed on the graph. 
Two classes of grid algorithms could be used. These are 
explained below. 

1 . The Wavefront Variant 

Starting with the source point at time zero, it 
progresses one step every time unit in all directions adding 
the appropriate nodes (these nodes which can be reached in the 
earliest at this time step) to a wavefront depending on the 
edge weights and direction of propagation. The wavefront keeps 
advancing every time unit until it hits the goal point. It is 
obvious that the progress in the inexpensive areas is bigger 
than expensive ones, since the length of step is proportional 
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Figure 7 : Neighbors In Grid Algorithm 
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to cost and time. Figure 8 shows an example in which the bigger 
strides take place in the low cost areas. 




One wavefront sweep will find the shortest path. But 
such shortest path may need a large number of time steps 
depending on edge weights. The number of nodes on the 
wavefront to be processed at every point in time varies; that 
causes difficulty in scheduling them on a distributed 
computer . 

2 . Relaxation-Based Approach 

The reason why this approach is called "relaxation- 
based" is because it is similar to Gauss-Seidel iteration 
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class of computation in Partial Differential Equations 
[HAB87] . 

Cij(m+ 1) = Cijim) +a)[c (l -_ 1) j(m) +c (i+1) j(m + 1) 

+ c ; , u _ 1} (m) + c ■ u+ 1} (m+ 1) -4c i j(m) ] 

where shows value of grid point in row i, column j of 

step m. 

The bracketed term indicates the change that occurs 
after each iteration as C£ f j(m) is updated to Ci,j (m+1) . 
Similarly, in our implementation, as shown in Figure 9, in 
each iteration, a node is updated. The cost of node is 
minimized by comparing it with the neighbors according to the 
following pseudocode: 

C current (4, j ) = min [ c old (i,j), c(i,j+l) + north_weight 

c(i+l,j) + east_weight 
c(i,j-l) + south_weight 
c(i-l,j) + west_weight ] 

We decided to implement the relaxation-based grid 
algorithm because of the straight-forward solutions on a 
distributed computer. 

a. Implementation 

The grid graph is represented as a two-dimensional 
array of records. Data structure and information carried by 
every node is shown in Figure 10 and the variables are 
described below. 
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Figure 9: A Node Described With Its Neighbors 
And Their Relations 



Current_cost : The beginning and updated cost value 
of the node. 

0ld_co8t : 01d_cost is the previous iteration cost 
value of the node and used for determining the change in cost. 
N,E,S,W: North, East, South and West neighbors. 
Weight: The weight on edge toward that direction. 
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Figure 10: Data Structure 

Distance : The step between two nodes (here it is 
always equal to one, but it changes in the modified 
algorithm) . 

Direction: String ( "North" , "East" , "South" , "West" ) 

and used for showing the path in the output . 

The pseudocode for the algorithm is presented in 
Figure 11. 

We obtained a data file which has 6400 raw data 
values representing the heights of 80*80 grid approximation 
of region. We constructed costs for the edges of the grids 
based on some function of these heights. For border nodes, 
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Algorithm: To find the Weighted Region Least Cost Path 



Given: Two dimensional grid of points: GRID 
Source point: (is, 3 s) 

Goal point: (ig, 3 g) 

Threshold: T 

Output : Cost of minimization from source to goal: Cost 
The minimum path from source to goal : 

(is, js) (ig, jg) 

procedure : 

Initialize 
Read the data file 
for all nodes loop 
calculate weight (i,j) 

end loop 

xx --At this point code will be inserted in the 

--modified version of the algorithm (see Fig. 14) 
flag = true 
while flag loop 
flag = false 
for all nodes loop 
find minimum value 

change = node ( i , j ) . old_cost-node ( i , j ) . current_cost 
if change > T then 
raise flag 

end if 
end loop 
end loop 

output the WRLCP 

Figure 11: The Psuedocode For Relaxation-Based Algorithm 



edges leading out of the grid where assigned very high 
positive values for weights, and minus one to distance. For 
all nodes we initialize the costs to a maximum positive number 
in order to use them in comparisons for finding the minimum. 

To output the Least Cost Path, we utilize a stack 
since the path traces back from the goal to the source point. 
So the output is displayed from source to goal as cost, 
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direction and distance (distance is equal to one in the 
straight-forward algorithm) are calculated at each step. 

A small-sized (4*4) example is presented in the 
Figure 12 for clear understanding. 

b. Parallelization 

In the relaxation-based approach, solutions can be 
straightforward on a vector computer or a on 2-d mesh of 
processors. But also there are some unattractive features such 
as the algorithm is intrinsically nonoptimal (that is true for 
the wavefront technique, too) . These algorithms do not reward 
the areas with the low cost by reducing the computations 
there. If the relaxation approach is used on a realistic 
number of processors, it is not known how to partition the 
computations. For instance, only one of the processors has the 
source point and starts computing, the computations of the 
rest of the processors for taking the minimum of infinite 
values are useless until they get the border values from that 
processor . 
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SOURCE 



Step 1. The heights are read from the 



file 



and loaded to nodes . 



Step 2 . The weights on the edges are calculated 
using some functions 
(e.g. abs ( 1024 - 1023) +1=2) 



Figure 12. a: Example For Relaxation-Based Algorithm 
The Values Inside The Nodes Are The 
Heights 
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Figure 12. b: Example For Relaxation-Based Algorithm 

The Values In The Nodes Are The Initial 
Costs From The Source. 
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GOAL 



12 . 



(continued from Figure 12. b) 

Step 4. This is the first iteration. In the 

latter iterations continue. 

Step 5. 14 becomes 12 after first iteration 

of cost minimization . 

Step 6. Thick arrows show the minimum least 
cost path 



Figure 12. c: Example For Relaxation-Based Algorithm 
The Values Inside The Nodes Are The 
Minimum Costs From The Source 
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III. A MODIFIED GRID ALGORITHM 



A. INTRODUCTION 

The complexity of the straight-forward grid algorithm in 

worst-case is 0(n 4 ) for an n x n grid. In the worst case, the 
WRLCP passes through all the nodes. Since the algorithm does 
not take the advantage of the low cost regions, we modified 
the algorithm to remedy this defficiency. 

B. VARIABLE GRID SIZE APPROACH 

We can use relaxation to obtain the shortest path from the 
source to different points without updating the nearest 
neihgbors as in the classical technique. We update the 
farthest neihgbors (two horizontal and two vertical in case 
of a 4-neihgbor technique) whose weighted distance is less 
than or equal to a prespecified constant. Determining such 
neighbors is done at the initialization phase. Rather than 
having the wavefront advance by a fixed distance, this 
algorithm has it advance by a fixed weight, thus bigger 
strides are made in cheaper regions. 

Information that decides which neighbor to propagate to is 
shown in Figure 13. It assumes that every grid point will have 
four records (N,E,S,W), and every record will have two 
components: the number of grid points to the neighbor to be 
updated (one in the classical method) , and the weighted 
distance to the neighbor. 
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Algorithm: I . To make horizantal search and elimination 
II. To decide which neghbor to propagate 

Given: Two dimensional grid of points: GRID (all edges 
Source point: (is,js) calculated) 

Goal point: (ig,jg) 

Stride 

Output : A searched grid and neighbor for node ( i , j ) to 
propogate 

procedure horizantal search and elimination 

for all rows i loop 

j = 1 

while j < jmax loop 

j = Do 

loop 

if weight from (is, j s ) to (ig,jg) < S then 

j = j+1 

end if 

until (weight from (is,js) to (ig,jg) => S) 
or (j = jmax) 

if weight from (is,js) to (ig,jg) > S then 

j = j-1 

end if 

delete all horizantal edges between (i,j 0 ) and 

( i , j ) 

( i , j o ) • east_weight = ( i , j ) . west_weight 

= weight difference 

( i , j o ) . east_distance = ( i , j ) . west_distance 

= j-Do 

end loop 
end loop 

procedure : proportion 

for all nodes (i,j) loop 
for north loop 

while { ( i , j ) . north_weight + ( i , j ) .north_weight 
( i , j ) . north_distance} < stride loop 
( i , j ) . north_distance : = ( i , j ) . north_distance + 1 
( i , j ) . north_weight : = old_value + 

(i, j+(i, j) . north_distance) . north_weight 

end loop 
end loop 
end loop 

repeat for other directions 



Figure 13: The Pseudocode for Modified Grid Algorithm 
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The procedure propagation assumes that the distance 
between successive grid points has been already computed; 
otherwise a more efficient strategy could be utilized. 

The modified relaxation algorithm is much faster than the 
traditional relaxation-based technique. The number of 
propagations in every polygon pj_ is reduced from Aj_ (the area 
of Pi) to Aj_/ (Wi**2) , where Wj_ is the weight of a unit distance 
in Pi. 

The user can change the threshold parameter in straight- 
forward algorithm, as he can also alter the stride parameter, 
in addition to the threshold parameter, in the modified 
algorithm. The threshold parameter ensures that the algorithm 
stops as it is with some threshold for the optimum. The stride 
parameter dictates the minimum weight of an edge thus making 
bigger strides in low cost areas. The stride parameter 
indirectly controls the number of edges and nodes to be 
eliminated, hence making the algorithm faster. 

An adventageous side-effect is the elimination of some 
nodes. These nodes are eliminated if they are not connected 
to more than one node (both source and goal points cannot be 
eliminated) . The psuedocode for these procedures are given 
in the same order as in the program: Horizantal search, edge 
eliminating, vertical search, edge eliminating and node 
eliminating in Figure 14. Furthermore, a continuation to the 
example in Figure 12 is given in Figure 15 . 
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Algorithm: To find the Weighted Region Least Cost Path 
in the modified version of relaxation-based 
grid algorithm 

Given: Two dimensional grid of points: GRID 
Source point: (is,js) 

Goal point: (ig,jg) 

Threshold: T 
Stride: S 

Output : Cost of minimization from source to goal: Cost 
The minimum path from source to goal: 

(is, js) (ig, jg) 

procedure : 

--This code is inserted at point xx in Figure 11 
Horizantal Search and Elimination (HSE) 

Vertical Search and Elimination ( similar to HSE) 
Eliminate the nodes which are not connected to more 
than one neighbor 

--The rest of the algorithm is the same except 
--while loop uses procedure propagation to propagate 



Figure 14 : The Pseudocode for Modified Grid Algorithm 



C. RESULTS 

We performed a series of experiments using the Meridian 
Ada compiler in Sparc stations in order to realize how the 
parameters effect time and cost. Some of results from these 
experiments are combined and shown in the Table I-IV. In Table 
I and II we tested threshold parameter and concluded that the 
greater the threshold value is, the faster the program is. 
Changing the threshold value does not have a significant 
effect on the least cost. The results taken from the 
experiments with modified program are presented in Table III- 
IV. In these experiments, we kept threshold value constant 
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(equal to zero) and changed the stride value between 1 and 
100. In the first row we realized that the cost of 
modifications increased the amount of time for computations. 
As seen form the last column, the time gets less even though 
the second trial took longer time. The number of the nodes 
increases depending on the stride value. Meanwhile, the 
disadvantage is that the least cost gets higher for higher 
stride values . 

D . DRAWBACKS 

The modified algorithm requires more overhead. The over- 
head as reflected by the results is worthwhile because the 
overall running time of the modified algorithm is smaller 
than the original algorithm for an approximately close 
stride value. Time overhead might be reduced by not requir- 
ing that distances between all neighboring gridpoint be pre- 
computed. Again, this algorithm is intrinsically suboptimal. 
There is a possibility that the graph is decomposed into more 
than one connected component due to the elimination of edges 
and nodes. If such is the case, then source and goal points 
might not be reachable from one another. This problem will 
manifest itself by having the minimum cost of the grid not 
changing from the assigned initial value. If there is more 
than one connected component, and the source and goal points 
are in the same connected component, then there should be a 
technique that avoids computations in the other connected 
components . 

The program codes for the modified grid algorithm are 
enclosed in Appendix B. 



25 




STEP 1. Horizontal search is performed 

(e.g. with stride = 10, in this case 
the weights an edges are added until 
it gets more than stride.) 

STEP 2. Edge Elimination is performed 



Figure 15. a: An Example For A Modified Algorithm (I) 
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GOAL 




SGSZaCB 

(continued from Figure 15. a) 

STEP 3 . Vertical Search is performed the same as 
horizontal search 

STEP 4. Edge eliminating is performed. 



Figure 15. b: An Example For A Modified Algorithm (II) 
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GOAL 




(continued from Figure 15. b) 

STEP 5. Node Elimination 

STEP 6. Cost Minimization and finding 

least cost path. 



Figure 15. c: An Example For A Modified Algorithm (III) 
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Table 1: EXPERIMENT SHOWING THE COST AND COMPUTATION TIME 
VERSUS THE THRESHOLD WITH CLASSICAL ALGORITHM FROM (1,1) TO 

(80,80) 



Threshold 


Cost 


Time(sec) 


0 


262 


6.9164 


10 


262 


5.3165 


20 


262 


2.8832 


30 


262 


2.8832 


50 


262 


2.8832 


100 


262 


2.8832 



Table 2: EXPERIMENT SHOWING THE COST AND COMPUTATION TIME 
VERSUS THE THRESHOLD WITH CLASSICAL ALGORITHM FROM (10,1) 

TO (55,80) 



Threshold 


Cost 


Time(sec) 


0 


308 


10.0996 


10 


308 


9.3163 


20 


308 


9.3163 


30 


308 


9.2996 


50 


308 


9.2996 


100 


308 


9.2330 
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Table 3: EXPERIMENT SHOWING THE COST AND COMPUTATION TIME 
VERSUS THE STRIDE WITH MODIFIED ALGORITHM FROM (1,1) TO (80,80) 



Threshold 


Stride 


Cost 


Time(sec) 


Number of 
Nodes 
Deleted 


0 


1 


262 


11.8329 


0 


0 


5 


276 


19.5826* 


1375 


0 


10 


298 


10.9996 


3177 


0 


20 


330 


9.3496 


4747 


0 


30 


356 


7.2497 


5391 


0 


50 


356 


4.0332 


5877 


0 


100 


356 


2.1666 


6172 



Table 4: EXPERIMENT SHOWING THE COST AND COMPUTATION TIME 
VERSUS THE STRIDE WITH MODIFIED ALGORITHM FROM (10,1) TO 

(55,80) 



Threshold 


Stride 


Cost 


Time(sec) 


Number of 
Nodes 
Deleted 


0 


1 


308 


17.3160 


0 


0 


5 


314 


19.6659* 


1375 


0 


10 


330 


10.9662 


3197 


0 


20 


336 


9.3996 


4730 


0 


30 


416 


7.797 


5390 


0 


50 


454 


4.0498 


5870 


0 


100 


454 


2.1499 


6168 



* : Overhead causes this increase in time 
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IV. THE ENVIRONMENT 



In this chapter we describe both the hardware and the soft- 
ware environments in which we implemented the parallel algo- 
rithm presented in Chapter V. The first section introduces the 
Transputer [TRA87], and the second section introduces Alsys 
Ada [ALS90] . A guide for program development using Ada on the 
Transputer is presented in the third section. 

A. TRANSPUTER 

The Transputer implementation is based on the concept of 
the Communicating Sequential Processes (CSP) . In order to uti- 
lize the transputer effectively, we need to understand the way 
it works. Transputer is a microprocessor with its own local 
memory storage and four links designed to communicate directly 
with other Transputers. The larger the number of processors 
in the network is, the more processing power, the more memory 
and links are available. The difficulties also grow with the 
number of processors. The most visible difficulty in the net- 
work is to avoid deadlock in which communication fails and 
results in processes waiting forever. 

There are different types of Transputer: T2 (T212, T222), 
T4 (T414, T425) AND T8 (T800, T801, T805) . We worked with T800 
Transputers. A block diagram of a T800 Transputer is presented 
in Figure 15. The major components of T800 Transputer, as 
seen, are memory, processor and communication system connect- 
ed via a 32 bit bus. 

The high level programming language OCCAM [OCC88] [OCC89] 
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is the primary language used for programming the Transputers. 
It is designed to run concurrent processes on a network of 
Transputers. Concurrence and communication are two main con- 
cepts in OCCAM. They allow processes to run simultaneously and 
transfer information through channels from process to pro- 
cess. Processes communicate by message passing, do not share 
variables, and synchronize only when they communicate. Commu- 
nication is synchronous and unbuffered. 

The host computer for the Transputer network that we use is a 
PC-286 . 
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B. ALSYS ADA 

Alsys Ada Compilation System, which consists of a compiler 
and a binder, is used for handling Ada programs in a 
Transputer programming environment. We used Alsys Ada as the 
language of choice for developing our WRLCP problem code. 
Below we describe the most salient feature that distinguishes 
Alsys Ada usage from ordinary programming in Ada. 

1. Channels 

Communication between Ada programs is provided by using 
transputer channels via the implementation defined package 
CHANNELS. The CHANNELS package contains a generic package 
CHANNEL_IO which defines input-output for values of a 
specified object type. READ, WRITE, READ_OR_FAIL , 
WRITE_OR_FAIL procedures within this generic package are used 
for input and output between channels. The distributed 
application is written as a set of independent programs for 
single or multiple Transputers and communicate through 
channels with unique names. 

We write a COMMON package, which contains declarations 
common to more than one Ada program and which also contains 
an instantiating of CHANNEL_IO to allow data to be 
communicated between independent programs . The data type of 
the channel, common to an application, is defined in a COMMON 
package. For example, to declare a channel of new NAMEl of 
type DATA_TYPE, we include the following in the COMMON 
package : 
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DATA_TYPE : (can be any type: integer, record...); 

package NAME1 is new CHANNELS . CHANNEL_IO (DATA_TYPE) ; 
In the programs it is used as follows: 
declaration-- Ada program segment begins here with declaration 
NAME 2 : CHANNELS .CHANNEL_REF : =CHANNELS . IN_PARAMETERS ( vi 
rtual channel number) ; 

--this channel is used for input from the other processors and virtual 
channel number is given by the programmer (e.g. 1 or 2 et:. 

NAME3 : CHANNELS . CHANNEL_REF : =CHANNELS . OUT_PARAMETERS (virtual 
channel number) ; --this is for output to the other processors 
begin --main program starts here 

NAME1 . READ (NAME2 , the same data_type) ;--get the input from the 
channel --NAME2 

NAME1 .WRITE (NAME3 , the same data_type) ; -- put the data_type 
value to --the channel NAME 3 . 

end 

2. Harnesses 

In order to run Ada programs in parallel on a single 
processor or a multi -transputer network, we need to use an 
interface, which is an occam process called a harness. A 
harness is used as a wrapping for the Ada program to be 
accepted by a Transputer. For every Ada program, two occam 
harnesses have to be created. Harnesses are explained in more 
details in section C. 

C. GUIDE FOR PROGRAM DEVELOPMENT 

In this section we present some helpful points to make 
program development easier. After learning the MAKE Program 
Maintenance Utility, the tools for checking the network, and 
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studying the examples, one will be ready for program 
development in this environment. 

1. MAKE program maintenance Utility 

MAKE is a utility program designed to help control of 
programming environment, to automate the process by 
determining which parts of the program is changed since the 
last compilation and rebuilds them accordingly. 

makefile is a script file, written by programmer and 
directs MAKE. You can find MAKE commands below: 

make family: Creates the Ada family and library sub- 
directories. This is a one-time-only operation. 

make: The standard command for building the executable 
programs after changes have been made to the source. 
make run: Executes the compiled program. 
make help: Displays the MAKE commands. 
make -n: Displays but do not execute commands. 
make check: Checks transputer topology. 
make clean: Deletes all files except source files. 
make *.o: Make Ada object codes. 

There is a batch file named doit_all that executes the 
first three commands -.make family, make, and make run. 

2 . Checking Tools 

Besides make check there are some more executable files 
which check the network topology: worm.exe, chknet.exe. 

3. Examples in Steps of Instructions 

It is better to start with complete examples to get used 

to it . 
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I. Make your own directory and copy the generic 
installation into it: 

>copy d:\al8ys037\8ource\generic\*.* . 

Complete documentation can be found in the files 

read. me and show. me. 

II. In this environment, there is an Ada source file 
proj.ada containing procedure proj. If you edit your own code 
with these names, it means you are ready to compile your code 
in Alsys Ada environment . 

III. You can type doit_all which executes three 
commands : 

make family, make, make run. 

IV. It is time to try the examples on a single 
transputer and then on multiple transputers. You can refer to 
the appendix C and Alsys Ada User Manuals. 

V. Communicating Ada processes on a single transputer 
needs the list of files below: 

makefile : A script file written by the programmer and 
executes the commands according to makefile. 

family. inv: This file creates the library environment 
(It does not change.). 

proj . inv: This file directs compiling and binding. 

main.occ : In order to integrate Ada with other 
languages a well defined interface is required. Ada programs 
may then be run in parallel on a single processor or 
distributed across a multi-transputer network, just as occam 
processes. This is a default occam harness provided as part 
of the compilation system in both source and compiled forms. 
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The main body of the harness consists of three processes 
operating in parallel: 

- A multiplexor which combines the error output and the 
standard output of the Ada program. 

- An error channel collector which collects any output 
from the error stream and routes it to the standard output 
stream of the server via the multiplexor. 

- A process which sets up the input and the output 
channel vectors of the Ada program and then invokes it, 
informing the other processes upon completion. 

merger. occ: This default harness is used to collect the 
error output from up to some number of Ada programs and send 
it to the standard output stream of the server (It does not 
change . ) . 

projh.occ : Each Ada (here PROJ.ADA) program has its own 
mini harness which provides a clean interface to the program 
in terms of the channels used. Main harness is used to invoke 
each of the mini harnesses in parallel. 

projh2 .occ : This is the dummy harness required to allow 
linking of a foreign Ada program with the occam libraries. 

main. Ink: Gives the file list to link. 

VI. The files needed for multiple transputers are 
mainly the same but they should be modified according to the 
network and presented below: 
makefile 

family. inv 

proj ? . inv: They should be as many as the number of 

transputers . 
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ma.inh.occ 



merger . occ 

proj?h.occ : They should be as many as the number of 
transputers . 

proj?h2.occ : They should be as many as the number of 

transputers . 

main .pgm: This is the only file not needed for one 

single transputer. It describes which virtual channels are 
equal to which physical channels. 

main . Ink 

An example is provided in Appendix C containing all 
these files, and it does not need to be changed for similar 
applications. In Figure 16 all these relations are shown. 
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common, ada 



proj 1 .ada 



projO.ada 



compile 


test_fam 


compile 




testjiD 


compile 









bind 



bine 



proj 1 .o 



projO.o 



projlh.occ occam 



projlh2.occ 



occart^ proj lh2. tax 



projlh.t8s 



ilink 



projlh.c8s 



projlh.mSs 



jlink 



iconf 



occam 



merger.occ 



merger.tSs 



projOh.occ occam 



projOh2.occ 



occarr^ projOh2.tax 



projOh.tSs _ 



occam 



mainh.occ 



mainh.tSs _ 



mainh.c8s 



mainh.m8s 



main.pgm 



lserver 



main.btl 



mainh.mSs 



Figure 17: Diagram Of The Steps Involved In Program Development 
Using Alsys Ada On Transputers 
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V. A MULTICOMPUTER ALGORITHM 



This chapter introduces a parallel algorithm for the WRLCP 
problem. In the first section we present the implementation 
on an INMOS Transputer network using Alsys Ada. In the second 
section we discuss a variant of that algorithm that uses less 
communication traffic. 

A. IMPLEMENTATION 

To simplify development of the algorithm, we made every 
possible effort to have it treat all processors symmetrically. 
This approach allows the algorithm to be scalable for a 
different number of processors without much change. We show 
the pseudocode for a version of the algorithm running on two 
processors in Figure 17 . 

The Root Transputer, which is the only Transputer having 
direct connection with the host PC, reads the data file and 
sends equal portions of the grid to other processors. At the 
beginning of the computation the Root Transputer sends the 
values of the threshold, the stride, the number of processors, 
the source and goal points to every processor. All the 
processors generate the weights on edges. All of them start 
cost minimization at the same time. Every processor updates 
the values of the grid points in its portion of the grid. For 
the points that are at the border of the grid portion assigned 
to a processor, the costs of the data obtained from the 
neighboring processor at the previous iteration are used. At 
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Algorithm: To find the Weighted Region Least Cost Path on 
two processors 

Given : Two dimensional grid of points: GRID 
Source point: (is, is) 

Goal point: (ig,jg) 

Threshold: T 

Output : Cost of minimization from source to goal: Cost 
The minimum path from source to goal : 

(is , js) (ig, jg) 

procedure: 

Initialize 

Re$d the dgta fiie 

write ( initial grid data) 

for all nodes loop 
calculate weight (i,j) 

end loop 

write (source, goal ,n,p, threshold) 
flag = true 
while flag loop 
flag = false 
for all nodes loop 

write (border costs) --exchanging the costs on border 
read (border costs) --nodes 
find minimum value 

change = node ( i , j ) . old_cost-node ( i , j ) . current_cost 
if change > T then 
raise flag 

end if 

end loop 

write (flag) --checking termination code 
read ( flag) 
end loop 

read (minimized costs) --for integration the solution, it 

--gets what the other processor did 

output the WRLCP 



Figure 18: Psuedocode For Parallel Algorithm 

From The One Processor's Point Of View. 

The Other Processor Will Have Read and Writes 
Interchanged . 



the end of the computation of an iteration, neighboring 
processors exchange values of the border grid points to be 
used in a later iteration. After every iteration, the 
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processors, as a whole check each other to decide whether to 
stop or to continue. If they decide that an adequate solution 
has been found they stop. All the processors write the 
obtained costs back to the Root Transputer which then displays 
the least cost path from the source to the goal points given. 

The pattern for "WRITE and READ" between the processors is 
shown in Figure 18. Other correct patterns exist, but it 
should be emphasized that a very important issue is deadlock 
avoidance . 



PROCESSOR- 1 


INITIAL GRID DATA 


PROCESSOR -2 






WRITE 

WRITE 




READ 

READ 


SOURCE, GOAL, N, P, THRESHOLD 










WRITE 

READ 




READ 
' WRITE 


EXCHANGING BORDER COSTS 






' FLAG ' 




WRITE 

READ 


READ 
' WRITE 


CHECKING TERMINATION CODE 
^ FLAG 






MINIMUM COST VALUES 




READ 




' WRITE 


INTEGRATING SOLUTION 



Figure 19: Write and Read Patterns 
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We have a running program on two processors, which is 
presented in appendix C. Furthermore, the methodology can be 
applied to more processors. Different patterns can be used to 
partition the grid graph on a number of processors . In Figure 
19 it is shown how the grid can be partitioned in different 
possible patterns, e.g. for four processors. 




B. A VARIANT ALGORITHM THAT SAVES ON COMMUNICATION 

The cost of communication can be very large since we have 
to exchange all the border values and checking parameters 
between the neighbor processors throughout every iteration of 
cost minimization. In reality, a lot of values exchanged 
across the border of the grid portions are redundant because 
they can remain unchanged through more than one iteration. 
That is why saving on communication becomes very important. 
As we try to speed up the algorithm, the cost of communication 
should not slow it down, especially if part of the 
communication is pragmatically useless. We present a modified 
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algorithm that takes advantage of the previous observation. 
In the modified approach, these border points which we changed 
in an iteration are marked. At the beginning of the data 
exchange between neighboring processors, each processor 
informs its neighbor about the number of the border points 
that changed. Then it proceeds to send the low cost values for 
only these points. More saving on communication can be 
achieved by using variant records. Unfortunately, we could not 
set this to work in the current development environment . The 
pseudocode for this approach is presented in Figure 20. 



Algorithm: The variant algorithm that saves on 
communication 

procedure ; 

for all border points node(i,j) loop 

change = node ( i , j ) . old_cost - node ( i , j ) . current_cost 
if change /= 0 then 
mark node ( i , j ) 

end if 
end loop 

count = number of marked nodes (i,j) 

write (count ) 

for 1 .. count loop 

write marked node(i,j) 

end loop 
read ( count ) 

for 1 .. count loop 

read marked node(i,j) 

end loop 



Figure 21: The Pseudocode For The Communication-Saver 
Variant Algorithm Of The Parallel Algorithm 
From The One Processor's Point Of View. 

The Other Processor Will Have Read and Writes 
Interchanged . 
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VI .CONCLUSIONS AND RECOMENDATIONS FOR FUTURE 

RESEARCH 



We have developed an efficient version of a grid based 
algorithm to solve the WRLCP problem. Our algorithm shows a 
significant decrease in the computation time in comparison to 
the original algorithm. The experienced loss in solution 
accuracy is not proportional to the saving in computation 
time . 

As a step towards developing a parallel version of this 
algorithm on a network of Transputers using Alsys Ada, we 
started developing simple parallel grid algorithm for the 
WRLCP problem in the above environment. Due to difficulties 
in dealing with that environment and due to lack of time, we 
stopped at the stage of developing adequate parallel 
algorithms for that environment, hoping that others will 
pursue our efforts towards the initial goal. 

The main emphasis in our parallel algorithm was compile 
time partitioning and mapping of data. Garcia [GAR89] 
implemented the parallel algorithm (Local, Asycnhronous and 
Iterative Parallel Procedures (LAIPP) Algorithm) presented in 
[SMI88] on a network of Transputers using Logical C. 
Scheduling was dynamic by farming out computations to 
available processors. Garcia's results showed that at a 
certain point, increasing the number of processors decreased 
the speedup. We attribute this to excessive communication 



45 



delays involved in the scheduling. We directed our effects 
towards nearest-neighbor patterns of communication, and we 
believe that this appropriate approach to handle this problem. 

Many problems need yet to be solved. As a starting point, 
software tools for automatic generation of Ada harnesses, and 
automatic mapping of Ada programs need to be acquired. This, 
and upgrade in the existing hardware setup will bring about 
more rapid program development. Currently, the process of 
program developing in that environment is extremely tedious. 

First order improvement to the existing parallel 
algorithms can be attained by using variant records for 
communication across channels. 

More serious improvement include: 

designing an asynchronous version of the parallel 
algorithm (data will be communicated only when needed) , 

- using queues or heaps can decrease the computations, 

- using the routing library developed in [FAL92] might 
provide us a way to compare our algorithm to more efficient 
algorithms that are not constrained to nearest neighbors. 

During the course of our work, we encountered problem with 
theoretical flavor which yet to be solved: 

- ensure that edge and node elimination in the algorithm 
(to reach a parallel analog of the modified sequential 
algorithm) does not separate the grid graph into different 
connected components . 
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APPENDIX A: SEQUENTIAL ALGORITHM SOURCE CODE 



--Title 

- -Author 

--Date 

--Revised 

--Course 

--Compiler 



STRAIGHT-FORWARD ALGORITHM 
CENGIZ EKIN 
20/06/92 
04/10/92 
THESIS 

MERIDIAN ADA 

--Description :Reads data from file, input starting and goal 
points --finds the minimum cost path, 
with TEXT_IO , OS_TYPES , TASK_CONTROL , CALENDAR; 
use TEXT_IO , OS_TYPES , CALENDAR; 



procedure MAIN1 is 

package INTEGER_INOUT is new INTEGER_IO ( INTEGER) ; 
package FLOAT_INOUT is new FLOAT_IO ( FLOAT) ; 
use FLOAT_INOUT, INTEGER_INOUT; 



START_TIME , 

END_TIME 

SX, SY, GX, GY, I, J 

Q 

VOLTA, VOLT 
NOP 

E , COUNTER 
El, SQ 

type ELER is 
record 
WEIGHT : INTEGER; 
DISTANCE : INTEGER; 
end record; 
type GRID_POINT is 
record 

CURRENT_COST , 

OLD_COST : INTEGER; 
N, 

E, 



: FLOAT ; 

; INTEGER ; 

: STRING ( 1. .1) := "y"; 

: INTEGER ; 

: INTEGER ; =80 ; 

: INTEGER; 

: STRING (1. .5) ; 
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s, 

W : ELER; 

DIRECTION : STRING ( 1 . . 5 ) ; 

end record; 

type GRID is array (0..(NOP +1 ) , 0 . . (NOP + 1 ) ) of GRID_POINT; 
B : GRID; 

INF : FILE_TYPE; 



type STORAGE is array (1..1000) of INTEGER; 
type ST0RAGE1 is array (1..1000) of STRING ( 1 .. 5 ) ; 
type STACK is 
record 



STORE 
LATEST 
end record; 
type STACK1 is 
record 
STORE 
LATEST 
end record; 



: STORAGE; 

: INTEGER := 0; 



: STORAGE 1 ; 

: INTEGER := 0; 



S : STACK; SI : STACK1 ; 



--This part is for the insertion of values into the STACK. 



procedure PUSH (S ; in out STACK; E : in INTEGER) is 
begin 



S. LATEST := S . LATEST + 1; 
S . STORE (S .LATEST) := E; 
end PUSH; 



procedure PUSHl (SI : in out STACK1 ; El : in STRING) is 



begin 



SI. LATEST := SI. LATEST + 1; 
SI . STORE (SI .LATEST) := El; 
end PUSHl; 



-- This part is to print the values from the STACK. 
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procedure POP (S : in out STACK; E : out INTEGER) is 
begin 

E : = S. STORE (S. LATEST ) ; 

S. LATEST := S . LATEST - 1; 
end POP; 



procedure POP1 (SI : in out STACKl ; El : out STRING) is 
begin 

El := SI. STORE (SI. LATEST) ; 

SI. LATEST := SI. LATEST - 1; 
end POP1 ; 



-- This function computes the execution time (CPU TIME) . 



function CLOK return FLOAT is 
function CLOCK return int ; 
pragma INTERFACE (C , CLOCK); 

T : int; 

S : FLOAT; 
begin 

task_control . pre_emption_of f ; 
T := CLOCK; 
if T = -1 then 

raise TIME_ERROR; 
else 

S := FLOAT (T) /I . 0E6; 
end i f ; 

task_control . pre_emption_on ; 
return S; 
end CLOK; 



procedure CAL_WEIGHT (I,J : in INTEGER) is 
begin 

If I = NOP then 

B ( I , J ) .N. WEIGHT := -1; 

B(I,J) .N. DISTANCE := -1; 

B ( I + 1 , J ) . CURRENT_COST :=10000; 
else 

B(I,J) .N. WEIGHT :=1+ ABS ( B ( I , J ) . CURRENT_COST - 
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B ( I + 1 , J ) . CURRENT_COST ) ; 
end i f ; 

If J = NOP then 

B ( I , J) .E. WEIGHT := -1; 

B(I,J) .E. DISTANCE := -1; 

B(I, J+l) . CURRENT_COST :=10000; 
else 

B(I, J) .E. WEIGHT : = 1 +ABS ( B ( I , J ) . CURRENT_COST - 
B ( I , J+l ) . CURRENT_COST) ; 
end i f ; 

If I = 1 then 

B ( I , J ) . S. WEIGHT : = -1; 

B(I,J) .S. DISTANCE := -1; 

B ( I - 1 , J ) . CURRENT_COST : =10000 ; 
else 

B(I, J) .S. WEIGHT : =1+ ABS ( B ( I , J ) . CURRENT_COST - B(I- 
1, J) . CURRENT_COST ) ; 
end i f ; 

If J = 1 then 

B ( I , J) .W. WEIGHT := -1; 

B(I, J) .W. DISTANCE := -1; 

B(I,J-1) . CURRENT_COST ;=10000; 
else 

B(I, J) .W. WEIGHT :=1+ ABS ( B ( I , J ) . CURRENT_COST - B(I,J- 
1) .CURRENT_COST) ; 
end i f ; 

end CAL_WEIGHT; 



procedure FIND_MIN (I,J : in INTEGER) is 
begin 

if B ( I , J ) . CURRENT_COST > ( B ( 1+ 1 , J ) . CURRENT_COST 



B ( I , J ) . CURRENT_COST 


. _ 


+abs (B (I, J) .N. WEIGHT) ) 

( B ( 1+1 , J) . CURRENT_COST 


then 


B (I, J) .DIRECTION 


. _ 


+abs (B ( I , J) .N. WEIGHT) ) 
" NORTH ■ ; 


/ 


end if; 

if B ( I , J ) . CURRENT_COST 


> 


( B ( I , J+l ) . CURRENT_COST 




B ( I , J ) . CURRENT_COST 


. _ 


+abs ( B ( I , J ) .E. WEIGHT) ) 

( B ( I , J+l ) . CURRENT_COST 


then 


B (I, J) .DIRECTION 


. 


+abs (B ( I , J) .E. WEIGHT) ) 
" EAST " ; 


/ 


end i f ; 

if B (I, J) .CURRENT_COST 


> 


( B ( I - 1 , J ) . CURRENT_COST 
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B ( I , J) 

B ( I , J) 
end if; 
if B ( I , J ) 

B ( I , J ) 

B ( I , J) 
end i f ; 
end FIND_MIN; 



. CURRENT_COST 
.DIRECTION 
. CURRENT_COST 
. CURRENT_COST 
.DIRECTION 



+abs ( B ( I , J) . S . WEIGHT) ) then 
( B ( I - 1 , J ) . CURRENT_COST 
+abs ( B ( I , J) . S. WEIGHT) ) ; 
"SOUTH" ; 

( B ( I , J-l ) . CURRENT_COST 
+abs (B(I,J) .W. WEIGHT) ) then 
( B ( I , J - 1 ) . CURRENT_COST 
+abs (B (I, J) .W. WEIGHT) ) ; 
"WEST 



-- Main program . . . 



begin 

while Q = "y" loop 
COUNTER := 1; 
for I in 1 . . NOP loop 



for J in 1 . . NOP 


loop 


B ( I , J ) .N. DISTANCE 


: =1 ; 


B ( I , J ) . E. DISTANCE 


: =1 ; 


B ( I , J ) . S. DISTANCE 


: =1 ; 


B ( I , J ) .W. DISTANCE 


: =1 ; 


end loop; 




loop; 





OPEN ( INF , MODE => IN_FILE , NAME => "ter.dat"); 

PUT ("THE DIMENSION OF MATRIX =" ); 

GET (NOP) ; 

PUT_LINE ("ENTER THE VOLTA (the optimization tolerance)"); 
PUT ("VOLTA = " ) ; GET (VOLTA); 

PUT_LINE ("ENTER THE SOURCE POINT !"); 

PUT ("SX = " ) ; GET (SX);PUT("SY = " );GET(SY); 

PUT_LINE ("ENTER THE GOAL POINT !"); 

PUT ( " GX = " ) ; GET (GX);PUT("GY = " ) ; GET (GY) ; 



-- CLOK function begins to compute the execution time. 
START_TIME := CLOK; 



--This part gets the heights from the file.. 



for ROW in 1..NOP loop 



54 



for COL in 1..N0P loop 

GET (INF, B (ROW, COL) . CURRENT_COST) ; 
end loop; 
end loop; 



--CLOK function finishes the computation of execution time. 
END_TIME := CLOK; 

PUT ("TOTAL TIME TO READ THE DATA FILE IS : "); 

PUT ( END_TIME-START_TIME ,4,4,0) ; new_line ; 

START_TIME := CLOK; 



--It determines the borders and calculates the weights of the 
edges . . 



for I in 1 . .NOP loop 
for J in 1..NOP loop 
C AL_WE I GHT ( I , J ) ; 
end loop; 
end loop; 



--It makes the costs max number in order to use them in 
comparisons for finding 
-- the minimum. . . . 



for I in 1 . . NOP loop 

for J in 1 . . NOP loop 

B (I, J) . CURRENT_COST := 10000; 

B ( I , J ) . OLD_COST : = B ( I , J ) . CURRENT_COST ; 
end loop; 
end loop ; 



--cost minimization. . . 



B(SX,SY) . CURRENT_COST := 0; 
while COUNTER > 0 loop 
COUNTER := 0; 
for I in 1..NOP loop 
for J in 1..NOP loop 
FIND_MIN ( I , J ) ; 

VOLT := B ( I , J ) . OLD_COST -B ( I , J ) . CURRENT_COST ; 
if VOLT > VOLTA then 

COUNTER := COUNTER +1; 
end i f ; 
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B ( I , J ) . OLD_COST : = B ( I , J ) . CURRENT_COST ; 
end loop; 
end loop; 
end loop; 



--output of least cost path.. 



loop 

PUSH (S, B (GX, GY) . CURRENT_COST) ; 
PUSHl (SI, B(GX, GY) .DIRECTION) ; 
SQ : =B (GX, GY) .DIRECTION; 
if SQ = "NORTH" then 

PUSH ( S , B (GX , GY ) .N. DISTANCE) ; 
GX := GX+1 ; 

elsif SQ = "EAST " then 

PUSH ( S , B (GX , GY) . E. DISTANCE) ; 
GY := GY + 1 ; 

elsif SQ = "SOUTH" then 

PUSH ( S , B (GX , GY) . S. DISTANCE) ; 
GX := GX-1 ; 

elsif SQ = "WEST " then 

PUSH (S, B (GX, GY) .W. DISTANCE) ; 
GY := GY - 1 ; 
else 
exit ; 
end if; 

exit when GX =SX and GY =SY ; 
end loop; 



DIRECTION" ) ; 



PUT_LINE ( " DISTANCE COST 

PUT_LINE (" 

loop 

PO P ( S , E ) ; PUT ( E ) ; PUT (" "); 

PO P ( S , E ) ; PUT ( E ) ; PUT (" "); 

POP1 (SI, El) ; 

if El = "NORTH" then 

PUT ( " SOUTH " ) ; new_l ine ; 
elsif El = "EAST " then 
PUT ( "WEST" ) ;new_line; 
elsif El = "SOUTH" then 
PUT ( "NORTH " ) ; new_line ; 
elsif El = "WEST " then 
PUT ( "EAST" ) ;new_line; 
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end i f ; 



exit when S . LATEST = 0; 
end loop; 



CLOSE (INF) ; 

--CLOK function finishes the computation of execution time. 
END_TIME := CLOK; 

PUT ("TOTAL TIME TO EXECUTE THE PROGRAM IS : "); 

PUT (END_.TIME-START_.TIME, 4,4,0) ;new_line; 

PUT ("RUN ONE MORE TIME :y(es) or n(o) :");GET(Q); 

if Q = "n" then exit; 

end if; 

end loop; 

end MAIN1; 
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APPENDIX B: MODIFIED ALGORITHM SOURCE CODE 



-Title 

-Author 

-Date 

-Revised 

-Course 

-Compiler 



MODIFIED ALGORITHM 

CENGIZ EKIN 

20/06/92 

04/10/92 

THESIS 

MERIDIAN ADA 



-Description rmodified program, 



searches and deletes nodes 



and edges 

with TEXT_IO,OS_TYPES, TASK_CONTROL , CALENDAR; 
use TEXT_IO, OS_TYPES, CALENDAR; 



procedure MAIN2 is 



package INTEGER_INOUT is new INTEGER_IO ( INTEGER) ; 
package FLOAT_INOUT is new FLOAT_IO ( FLOAT) ; 
use FLOAT_INOUT, INTEGER_INOUT; 



START_TIME , 

END_TIME : FLOAT; 



. . 5 ) ; 

. .1) := "y" ; 

: = 8 0 ; - - 5 0 0 ; 

type ELER is 
record 

WEIGHT : INTEGER; 

DISTANCE : INTEGER; 
end record; 
type GRID_POINT is 
record 

CURRENT_COST , 

OLD_COST : INTEGER; 



SX, SY,GX ( GY, Z , C02 , C03 , C04 , C05 
E, COUNTER, I, J,L 
El, SQ 
Q 

DELI, A, VOLTA, VOLT 

NOP 

MARK 



INTEGER; 

INTEGER; 

STRING ( 1 

STRING ( 1 

INTEGER; 

INTEGER 

BOOLEAN; 
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N, 

E, 

S, 

W : ELER ; 

DIRECTION : STRING ( 1 .. 5 ) ; 

ACTIVE : BOOLEAN; 

end record; 

type GRID is array (0..(NOP +1 ) , 0 . . (NOP+1 ) ) of GRID_POINT; 
B .-GRID; 

INF : FILEJTYPE; 



type STORAGE is array (1..1000) of INTEGER; 
type S TO RAGE 1 is array (1..1000) of STRING ( 1 . 
type STACK is 
record 



STORE 
LATEST 
end record; 
type STACK1 is 
record 
STORE 
LATEST 
end record; 



: STORAGE; 

: INTEGER := 0; 



: STORAGE 1 ; 

: INTEGER := 0; 



5) ; 



S : STACK; SI : STACK1 ; 



--This part is for the insertion of values into the STACK, 
procedure PUSH (S : in out STACK; E : in INTEGER) is 
begin 

S. LATEST := S . LATEST + 1; 

S . STORE (S .LATEST) := E; 
end PUSH; 



procedure PUSH1 (SI : in out STACKl ; El : in STRING) is 
begin 



SI. LATEST := Si. LATEST + 1; 
SI. STORE (SI. LATEST) := El; 
end PUSH1 ; 
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-- This part is to print the values from the STACK. 



procedure POP (S : in out STACK; E : out INTEGER) is 
begin 

E := S. STORE (S. LATEST ) ; 

S. LATEST := S. LATEST - 1; 
end POP; 



procedure P0P1 (SI : in out STACK1; El : out STRING) is 
begin 

El := SI .STORE (SI .LATEST) ; 

SI. LATEST := SI. LATEST - 1; 
end POP1; 



-- This function computes the execution time (CPU TIME) . 



function CLOK return FLOAT is 
function CLOCK return int ; 
pragma INTERFACE (C, CLOCK); 

T : int; 

S : FLOAT; 

begin 

task_control . pre_empt ion_of f ; 

T := CLOCK; 

if T = -1 then 

raise TIME_ERROR; 
else 

S := FLOAT (T) /I .0E6; 
end i f ; 

task_control .pre_emption_on; 
return S; 
end CLOK; 



procedure CAL_WEIGHT (I,J : in INTEGER) is 
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begin 

If I = NOP then 

B ( I , J ) .N. WEIGHT := -1; 

B(I, J) .N. DISTANCE := -1; 

B(I+1, J) . CURRENT_COST :=10000; 
else 

B(I,J) .N. WEIGHT : =1+ ABS ( B ( I , J ) . CURRENT_COST - 
B ( 1+1 , J) . CURRENT_COST) ; 
end i f ; 

If J - NOP then 

B(I,J) .E. WEIGHT := -1; 

B(I,J) .E. DISTANCE := -1; 

B(I,J+1) . CURRENT_COST :=10000; 
else 

B (I, J) .E. WEIGHT := 1+ABS (B ( I , J) . CURRENT_COST - 
B ( I , J+l ) . CURRENT_COST) ; 
end if; 

If I = 1 then 

B ( I , J) .S. WEIGHT := -1; 

B(I,J) .S. DISTANCE := -1; 

B ( I - 1 , J ) . CURRENT_COST :=10000; 
else 

B ( I , J ) . S . WEIGHT :=1+ ABS { B ( I , J ) . CURRENT_COST - 
1, J) . CURRENT_COST ) ; 
end if; 

If J = 1 then 

B ( I , J ) .W. WEIGHT := -1; 

B(I, J) .W. DISTANCE := -1; 

B(I,J-1) . CURRENT_COST ;=10000; 
else 

B ( I , J ) .W. WEIGHT :=1+ ABS(B(I,J) . CURRENT_COST - E 
1) . CURRENT_COST ) ; 
end if; 

end CAL_WEIGHT; 



procedure FIND_MIN (I,J : in INTEGER) is 
begin 

if B ( I , J ) .ACTIVE = TRUE then 
if B ( I , J ) .N. DISTANCE /= -1 then 
if B ( I , J ) . CURRENT_COST > 

( B ( ( I+abs ( B ( I , J ) .N. DISTANCE) ) ,J) . CURRENT_COST 

+abs (B ( I , J) .N. WEIGHT) ) 

B ( I , J ) . CURRENT_COST : = 



B ( I - 



(I, J- 



then 
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( B ( (I+abs (B(I, J) .N. DISTANCE) ) , J) . CURRENT_COST 

+abs ( B ( I , J ) .N. WEIGHT) ) ; 
B(I,J) .DIRECTION := "NORTH"; 

end if; 
end i f ; 

if B (I , J) .E. DISTANCE /= -1 then 
if B(I,J) . CURRENT_COST > 

( B ( I , (J+abs ( B ( I , J) .E. DISTANCE) ) ) . CURRENT_COST 

+abs ( B ( I , J ) . E . WEIGHT) ) then 

B ( I , J ) . CURRENT_COST : = 

( B ( I , (J+abs ( B ( I , J) .E. DISTANCE) ) ) . CURRENT_COST 

+abs ( B ( I , J ) .E. WEIGHT) ) ; 
B(I,J) .DIRECTION : = "EAST "; 

end i f ; 
end if; 

if B ( I, J) .S. DISTANCE /= -1 then 
if B ( I , J ) . CURRENT_COST > (B((I- 

abs ( B ( I , J) . S. DISTANCE) ) , J) . CURRENT_COST 

+abs ( B ( I , J) . S .WEIGHT) ) then 
B ( I , J ) . CURRENT_COST := (B((I- 
abs ( B ( I , J ) . S. DISTANCE) ) , J) . CURRENT_COST 

+abs ( B ( I , J) . S. WEIGHT) ) ; 
B(I,J) .DIRECTION := "SOUTH"; 

end i f ; 
end i f ; 

if B ( I , J ) .W. DISTANCE /= -1 then 
if B ( I , J ) . CURRENT_COST > (B(I,(J- 

abs ( B ( I , J) .W. DISTANCE) ) ) . CURRENT_COST 

+abs(B(I, J) .W. WEIGHT) ) then 
B ( I , J ) . CURRENT_COST := (B(I,(J- 
abs ( B ( I , J) .W. DISTANCE) ) ) . CURRENT_COST 

+abs ( B ( I , J) .W. WEIGHT) ) ; 

B ( I, J) .DIRECTION := "WEST "; 

end i f ; 
end i f ; 
end i f ; 
end FIND_MIN; 



-- Main program . . . 



begin 

while Q = "y" loop 
MARK := TRUE; 
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COUNTER := 1; 



for I in 1 
for J in 
B ( I , J ) 
B ( I , J ) 
B ( I , J ) 
B ( I , J ) 
B ( I , J) 
end loop; 
end loop; 



. NOP loop 
1 . . NOP loop 



.N. DISTANCE 
.E. DISTANCE 
. S. DISTANCE 
.W. DISTANCE 
. ACTIVE 



= 1 ; 

= 1 ; 

= 1 ; 

= 1 ; 

: =TRUE ; 



OPEN ( INF , MODE => IN_FILE, NAME => "ter.dat"); 

PUT ("THE DIMENSION OF MATRIX =" ); 

GET (NOP) ; 

PUT_LINE ("ENTER THE DELTA (the jumping value)."); 

PUT ("DELI = " ) ; GET (DELI); 

PUT_LINE ("ENTER THE VOLTA (the optimization tolerance) .") ; 
PUT ("VOLTA = " ) ; GET (VOLTA); 

PUT_LINE ("ENTER THE SOURCE POINT !"); 

PUT ("SX = " ) ; GET (SX);PUT("SY = ");GET(SY); 

PUT_LINE ("ENTER THE GOAL POINT !"); 

PUT ( " GX = " ) ; GET (GX);PUT("GY = " ) ; GET ( GY) ; 



-- CLOK function begins to compute the execution time. 
START_TIME : = CLOK; 



--This part gets the heights from the file.. 



for ROW in 1..NOP loop 
for COL in 1..NOP loop 

GET (INF, B (ROW, COL) . CURRENT_COST) ; 
end loop; 
end loop; 



--CLOK function finishes the computation of execution time. 
END_TIME : = CLOK ; 

PUT ("TOTAL TIME TO READ THE DATA FILE IS : "); 

PUT ( E ND_T I ME - S T ART_T I ME ,4,4,0); new_line ; 

START_TIME := CLOK; 



--It determines the borders and calculates the weights of the 
edges . . 
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for I in 1..N0P loop 

for J in 1..N0P loop 
CAL_WEIGHT ( I , J) ; 
end loop; 
end loop; 



--It makes the costs max number in order to use them in 
comparisons --for finding the minimum. . . . 



for I in 1 . . NOP loop 

for J in 1 . . NOP loop 

B ( I , J ) . CURRENT_COST : = 10000; 

B ( I , J ) . OLD_COST : = B ( I , J ) . CURRENT_COST ; 
end loop; 
end loop; 



--horizontal search.. 



B(SX,SY) . CURRENT_COST := 0; I:=1;J:=1; 
loop 

A : = B (I, J) .E. WEIGHT; 

C02 : = 0 ; 

if A <= DELI then 
loop 

exit when (A > DELI) ; 
if ( J>=NOP) then exit; end if; 
if ( ( I=SX and J=SY) or (I=GX and J=GY) ) and (MARK =TRUE ) 

then 

MARK := FALSE; exit; end if; 

MARK := TRUE; 

J := J + 1; 

C02 := C02 + 1; 

A : = A + B(I, J) .E. WEIGHT; 
end loop; 
if C02 > 0 then 

B(I,J-C02) .E. WEIGHT := A - B ( I , J ). E . WEIGHT ; 

B ( I, J-C02 ) .E. DISTANCE := C02 ; 

B(I, J) .W. WEIGHT := B(I,J-C02).E.WEIGHT; 

B ( I, J) .W. DISTANCE := B ( I , J-C02 ). E . DISTANCE ; 

end i f ; 
else 

J : = J + 1; 
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end i f ; 

if J - NOP then 
I := I + 1; 

J := 1; 
end 1 f ; 

exit when I = NOP + 1; 
end loop; 



-- horizontal edge eliminating... 



C03 := 0; J := 1;I ;=1; 
loop 

if B ( I, J) .E. DISTANCE > 1 then 
Z := B ( I, J) .E. DISTANCE + J; 
loop 

J := J+l ; 

exit when J = Z; 

B ( I, J) .E. DISTANCE := -1; 
B(I, J) .W. DISTANCE := -1; 
C03 := C03 + 2; 
end loop ; 
else 

J := J+l; 
end i f ; 

if J = NOP then 
I := I + 1; 

J := 1; 
end i f ; 

exit when I = NOP + 1; 
end loop; 



--vertical search. . . 



I := 1 ; J : = 1 ; MARK TRUE; 
loop 

A := B ( I , J) . N. WEIGHT; 

C02 : = 0 ; 

if A <= DELI then 
loop 

exit when (A > DELI) ; 
if (I>=NOP) then exit; end if; 
if ( ( I=SX and J=SY) or (I-GX and J=GY) ) and (MARK =TRUE) 

then 
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MARK := FALSE; exit; end if; 
MARK := TRUE; 

I := I + 1; 

C02 ;= C02 + 1; 

A : = A + B(I, J) .N. WEIGHT; 



end loop ; 
if C02 > 0 then 
B ( I -C02 , J) .N. WEIGHT 
B ( I-C02 , J) .N. DISTANCE 
B ( I , J ) . S. WEIGHT 
B ( I, J) . S .DISTANCE 
end i f ; 
else 



A - B ( I , J) .N. WEIGHT; 
C02 ; 

B ( I-C02 , J) . N .WEIGHT; 

B ( I-C02 , J) .N. DISTANCE; 



I := I + 1; 
end i f ; 

if I = NOP then 
J ; = J + 1 ; 

I := 1; 
end i f ; 

exit when J = NOP + 1; 
end loop; 



--vertical edge eliminating... 



I := 1 ; J : = 1 ; 
loop 

if B ( I, J) .N. DISTANCE > 1 then 
Z ;= B ( I , J ) .N. DISTANCE + I; 
loop 

I := 1+1 ; 

exit when I=Z; 

B(I, J) .N. DISTANCE := -1; 
B ( I, J) .S. DISTANCE := -1; 
C03 ;= C03 + 2; 
end loop; 
else 

I := 1+1; 
end i f ; 

if I = NOP then 
J := J + 1; 

I : = 1 ; 
end i f ; 

exit when J = NOP + 1; 
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end loop; 



--node eliminating . . . 



C05 := 0 ; J :=1;I:=1; 
for I in 1 . . NOP loop 
if not ( I=SX and l=SY) and 
C04 : = 0 ; 

if B(I, I) .N. DISTANCE 
C04 := C04 + 1; 

end if; 

if B(I, I) . E. DISTANCE 
C04 := C04 + 1; 
end i f ; 

if B ( I , I ) . S. DISTANCE 
C04 := C04 + 1; 
end if; 

if B ( I , I ) .W. DISTANCE 
C04 := C04 + 1; 
end if; 

if C04 >= 3 then 



not(I=GX and I 

< 1 then 

< 1 then 

< 1 then 

< 1 then 



B (I, I) .ACTIVE := FALSE; 
if Bd, I) .N. DISTANCE /= -1 then 
B ( ( I+abs ( B ( I , I ) .N. DISTANCE) ) , I 
B(I, I) .N. DISTANCE := -1; 

C03 : =C03 + 1; 
end i f ; 



if B (I, I) .E. DISTANCE /= -1 then 
B ( I , (I+abs (B(I, I) .E. DISTANCE) ) 
B(I, I) .E. DISTANCE ;= -1; 

C03 : =C03 + 1; 
end if; 

if B (I, I) .S. DISTANCE /= -1 then 
B ( ( I-abs (B ( I , J) . S. DISTANCE ) ) ,J 
B(I, I) .S. DISTANCE := -1; 

C03 : =C03 + 1; 
end if; 

if B(I, I) .W. DISTANCE /= -1 then 
B ( I , (I-abs (B(I, I) .W. DISTANCE) ) ) 
B(I, I) .W. DISTANCE ;= -1; 

C03 : =C03 + 1; 
end i f ; 



=GY) then 



) . S .DISTANCE : = -1 ; 



) .W. DISTANCE ;= -1; 



) .N. DISTANCE : = -1; 



. E. DISTANCE := -1; 
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C05 := C05 +1; 
end i f ; 
end i f ; 

for J in 1+1 . . NOP loop 
C04 : = 0 ; 

if not(I=SX and J=SY) and 
if B(I, J) .N. DISTANCE < 1 
C04 := C04 + 1; 
end i f ; 

if B ( I , J ) .E. DISTANCE < 1 
C04 := C04 + 1; 
end i f ; 

if B(I, J) .S. DISTANCE < 1 
C04 := C04 + 1; 
end i f ; 

if B ( I , J ) .W. DISTANCE < 1 
C04 := C04 + 1; 
end i f ; 

if C04 >= 3 then 



not (I=GX 
then 



then 



then 



then 



and J=GY ) 



then 



B ( I , J ) .ACTIVE := FALSE; 

if B ( I , J ) .N. DISTANCE /= -1 then 

B ( (I+abs (B(I, J) .N. DISTANCE) ) , J) .W. DISTANCE 
B (I, J) .N. DISTANCE := -1; 

C03 :=C03 + 1; 
end i f ; 

if B ( I , J ) .E. DISTANCE /= -1 then 

B ( I , (J+abs (B(I, J) . E. DISTANCE) ) ) .W. DISTANCE 
B ( I, J) .E. DISTANCE := -1; 

C03 : =C03 + 1; 
end i f ; 

if B(I,J) .S. DISTANCE /= -1 then 

B ( (I-abs (B(I, J) . S. DISTANCE) ),J) ,N. DISTANCE 
B (I, J) .S. DISTANCE := -1; 

C03 : =C03 + 1; 
end i f ; 

if B(I, J) .W. DISTANCE /= -1 then 

B(I, (J-abs (B(I, J) ,W. DISTANCE) ) ) .E. DISTANCE 
B(I, J) .W. DISTANCE := -1; 

C03 :=C03 + 1; 
end i f ; 

COS := C05 +1; 
end i f ; 
end i f ; 
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C04 := 0; 

if not ( J = SX and I = SY) 
if B(J, I) .N. DISTANCE 
C04 : = C04 + 1 ; 
end i f ; 

if B ( J , I ) . E. DISTANCE 
C04 := C04 + 1; 
end i f ; 

if B ( J, I) .S. DISTANCE 
C04 := C04 + 1; 
end i f ; 

if B(J, I) .W. DISTANCE 
C04 := C04 + 1; 

end i f ; 

if C04 >= 3 then 



and not(J=GX and I=GY) 
< 1 then 



< 1 then 



< 1 then 



< 1 then 



then 



B(J, I) .ACTIVE := FALSE; 
if B(J, I) .N. DISTANCE /= -1 then 
B ( ( J+abs (B(J, I) .N. DISTANCE) ) , I) . S .DISTANCE : 
B(J, I) .N. DISTANCE := -1; 

C03 : =C03 + 1; 
end i f ; 

if B(J, I) .E. DISTANCE /= -1 then 

B ( J, ( I+abs ( B ( J, I) . E. DISTANCE) ) ) .W. DISTANCE 
B(J, I) .E. DISTANCE := -1; 

C03 : =C03 + 1; 
end i f ; 

if B(J,I) .S. DISTANCE /= -1 then 

B ( (J-abs (B(J, I) . S. DISTANCE) ) ,1) .N. DISTANCE 
B(J,I) .S. DISTANCE := -1; 

C03 : =C03 + 1; 
end i f ; 

if B(J, I) .W. DISTANCE /= -1 then 

B ( J , ( I-abs ( B ( J, I) .W. DISTANCE) ) ) . E. DISTANCE 
B (J, I) .W. DISTANCE := -1; 

C03 : =C03 + 1; 
end i f ; 



C05 := C05 +1; 
end i f ; 
end if; 
end loop; 
end loop; 
put (co5) ; 



- 1 ; 



- 1 ; 



- 1 ; 
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--cost minimization... 



while COUNTER > 0 loop 
COUNTER := 0; 
for I in 1..NOP loop 
for J in 1..NOP loop 
FIND_MIN ( I , J ) ; 

VOLT := B ( I , J ) . OLD_COST - B ( I , J ) . CURRENT_COST ; 
if VOLT > VOLTA then 

COUNTER := COUNTER +1; 
else 

COUNTER := COUNTER; 
end i f ; 

B ( I , J ) . OLD_COST : = B ( I , J ) . CURRENT_COST ; 
end loop; 
end loop; 
end loop; 



--output of least cost path. . 



loop 

PUSH ( S , B ( GX , GY ) . CURRENT_COST ) ; 
PUSH1 (SI, B(GX, GY) .DIRECTION) ; 

SQ : =B (GX, GY) .DIRECTION; 
if SQ = "NORTH" then 

PUSH ( S , B ( GX , GY ) .N. DISTANCE) ; 
GX := GX+B (GX, GY) .N. DISTANCE; 
elsif SQ = "EAST * then 

PUSH ( S , B ( GX , GY ) . E. DISTANCE) ; 
GY := GY+B ( GX, GY ) .E. DISTANCE; 
elsif SQ = "SOUTH" then 

PUSH ( S , B ( GX , GY ) .S. DISTANCE) ; 
GX := GX-B(GX, GY) .S. DISTANCE; 
elsif SQ = "WEST “ then 

PUSH (S, B(GX,GY) .W. DISTANCE) ; 
GY := GY-B(GX, GY) .W. DISTANCE; 
else 
exit ; 
end i f ; 

exit when GX =SX and GY =SY ; 
end loop; 



PUT_LINE ( " DISTANCE COST DIRECTION"); 
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PUT_LINE(" "); 

loop 

POP ( S , E ) ; PUT ( E) ; PUT ( " " ) ; 

POP ( S , E ) ; PUT ( E ) ; PUT (" "); 

POPl (SI, El) ; 

if El = "NORTH" then 

PUT ( " SOUTH " ) ; new_l ine ; 
elsif El = "EAST " then 
PUT ( "WEST" ) ;new_line; 
elsif El = "SOUTH" then 
PUT ( "NORTH " ) ; new_line ; 
elsif El = "WEST " then 
PUT ( " EAST " ) ; new_l ine ; 
end i f ; 

exit when S. LATEST =0; 
end loop; 



CLOSE (INF) ; 

--CLOK function finishes the computation of execution time. 
END_TIME := CLOK; 

PUT ("TOTAL TIME TO EXECUTE THE PROGRAM IS : "); 

PUT ( END_TIME-START_TIME , 4,4, 0) ;new_line; 

PUT ("RUN ONE MORE TIME :y(es) or n(o) :");GET(Q); 

if Q = "n" then exit; 

end i f ; 

end loop; 

end MAIN2; 
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APPENDIX C: PARALLEL ALGORITHM SOURCE CODE 



--Title 
- -Author 
- -Date 
- -Revised 
--Course 
- -Compiler 



PRO JO .ADA 
CENGIZ EKIN 
20/06/92 
04/10/92 
THESIS 
ALSYS ADA 

--Description : Reads data from file, input 
-- points, finds the minimum cost path. 



starting and goal 



with TEXT_IO , COMMON, CHANNELS; 
use COMMON; 
procedure PROJO is 

package INTEGER_INOUT is new TEXT_IO . INTEGER_IO ( INTEGER) ; 
package FLOAT_INOUT is new TEXT_IO . FLOAT_IO ( FLOAT ) ; 
use FLOAT_INOUT, INTEGER_INOUT ; 

INF : TEXT_IO.FILE_TYPE; 

B : GRID ; GRI ;GRID_POINT; 



counter : integer :=1; 

-- communication channels that are used 



OutToMars : CHANNELS . CHANNEL_REF := 
CHANNELS . OUT_PARAMETERS ( 2 ) ; 

InFromMars : CHANNELS . CHANNEL_REF : = 
CHANNELS . IN_PARAMETERS ( 2 ) ; 



--This procedure calculates the weights on edges and builds 
up an hyphotetical 

--wall for processes to stay in the terrain... 



procedure CAL_WEIGHT (I,J : in INTEGER) is 
begin 

If I = 10 then 

B ( I , J) .N. WEIGHT : = -1; 

B(I, J) .N. DISTANCE := -1; 

B ( I + 1 , J ) . CURRENT_COST .-=10000 ; 
else 

B (I, J) .N. WEIGHT : =1+ ABS ( B ( I , J ) . CURRENT_COST - 
B ( 1+1 , J) . CURRENT_COST ) ; 
end if; 

If J = 10 then 
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B ( I , J ) .E. WEIGHT := -1; 

B(I,J) .E. DISTANCE := -1; 

B(I, J+l) .CURRENT_COST :=10000; 
else 

B(I,J) .E. WEIGHT : = 1+ABS ( B ( I , J) . CURRENT_COST - 
B(I, J+l) . CURRENT_COST) ; 
end i f ; 

If I - 1 then 

B ( I , J) . S. WEIGHT := -1; 

B(I,J) .S. DISTANCE := -1; 

B ( I - 1 , J ) . CURRENT_COST :=10000; 
else 

B(I, J) .S. WEIGHT : =1+ ABS ( B ( I , J ) . CURRENT_COST - 
B ( 1-1 , J) . CURRENT_COST ) ; 
end i f ; 

If J = 1 then 

B(I,J) .W. WEIGHT := -1; 

B(I, J) .W. DISTANCE := -1; 

B(I, J-l) . CURRENT_COST :=10000; 
else 

B(I, J) .W. WEIGHT :=1+ ABS ( B ( I , J ) . CURRENT_COST - 
B(I, J-l) . CURRENT_COST) ; 
end i f ; 

end CAL_WEIGHT; 



-- This part finds the minimum cost between current node and 
its four 

--neighbors (north, east , south, west ) 



procedure FIND_MIN (I,J : in INTEGER) is 
begin 



if B ( I , J ) . CURRENT_COST > 

B ( I , J ) . CURRENT_COST : = 

B ( I , J ) .DIRECTION := 

end if; 

if B ( I , J ) . CURRENT_COST > 

B ( I , J ) . CURRENT_COST : = 

B ( I, J) .DIRECTION := 

end i f ; 



( B ( 1+1 , J) . CURRENT_COST 
+abs (B(I,J) .N. WEIGHT) ) then 
( B ( 1+1 , J ) . CURRENT_COST 
+abs (B (I, J) .N. WEIGHT) ) ; 

" NORTH " ; 

( B ( I , J+l ) . CURRENT_COST 
+abs (B (I, J) .E. WEIGHT) ) then 
(B ( I , J+l ) . CURRENT_COST 
+abs (B (I, J) .E. WEIGHT) ) ; 

" EAST " ; 
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if B (I, J) ,CURRENT_COST > 

B ( I , J ) . CURRENT_COST : = 

B ( I, J) .DIRECTION := 

end i f ; 

if B ( I , J ) . CURRENT_COST > 

B(I,J) . CURRENT_COST : = 

B(I,J) .DIRECTION := 

end i f ; 
end FIND_MIN; 



( B ( I - 1 , J ) . CURRENT_COST 
+abs (B(I,J) .S. WEIGHT) ) then 
( B ( I - 1 , J ) . CURRENT_COST 
+abs ( B ( I , J ) . S. WEIGHT) ) ; 
"SOUTH" ; 

( B ( I , J-l ) . CURRENT_COST 
+abs (B(I,J) .W. WEIGHT) ) then 
(B(I, J-l) . CURRENT_COST 
+abs ( B ( I , J ) .W. WEIGHT) ) ; 
"WEST " ; 



begin 

TEXT_IO . OPEN (INF, MODE => TEXT_IO . IN_FILE, NAME => 

" ter . dat " ) ; 

TEXT_IO . PUT_LINE ( " THE DIMENSION OF MATRIX ="); 
INTEGER_INOUT . GET (N) ; 

TEXT_IO . PUT_LINE ( " THE NUMBER OF PROCESSORS = "); 
INTEGER_INOUT . GET ( P ) ; 

TEXT_IO . PUT_LINE ("ENTER THE OPTIMIZATION TOLERANCE"); 
INTEGER_INOUT . GET (VOLTA) ; 

TEXT_IO.PUT_LINE ("ENTER THE SOURCE POINT !"); 

TEXT_IO . PUT_LINE ( " SX = " ) ; INTEGER_INOUT . GET (SX); 
TEXT_IO . PUT_LINE ( " SY = " ) ; INTEGER_INOUT . GET ( SY) ; 
TEXT_IO . PUT_LINE ("ENTER THE GOAL POINT !"); 

TEXT_IO . PUT_LINE ( "GX = " ) ; INTEGER_INOUT . GET (GX); 
TEXT_IO. PUT_LINE ( "GY = ") ; INTEGER_INOUT . GET ( GY) ; 



--This part gets the heights from the file.. 



for ROW in 1..10 loop 
for COL in 1..10 loop 

INTEGER_INOUT . GET (INF, B ( ROW, COL ) . CURRENT_COST) ; 
end loop; 
end loop; 



--It passes the heights to the other processors.. 



for I in 1 . . 5 loop 
for J in 1 . .10 loop 
GRI := B ( I , J ) ; 



74 



DATA_IO. WRITE (OutToMars , GRI ) ; 
end loop; 
end loop; 



--It determines the borders and calculates the weights of the 
edges . . 



for I in 6 . . 10 loop 
for J in 1. .10 loop 
CAL_WEIGHT(I, J) ; 
end loop; 
end loop; 



--It makes the costs max number in order to use them in 
comparisons for finding 
-- the minimum. . . . 



for I in 6 . . 10 loop 

for J in 1 . . 10 loop 

B ( I , J ) . CURRENT_COST := 10000; 

B ( I , J ) . OLD_COST := B ( I , J ) . CURRENT_COST ; 
end loop; 
end loop; 



--This part sends dim of matrix, no of 
proccessors , volta , sourceand goal points. 



GRI . CURRENT_COST : =N; GRI . OLD_COST : =P ; GRI . N . WEIGHT : =VOLTA; 

GRI . E . WEIGHT ; =SX ; GRI . S . WEIGHT : =SY ; GRI . W . WEIGHT : =GX ; GRI . E . DIS 
TANCE : =GY; 

DATA_IO. WRITE (OutToMars, GRI) ; 



--cost minimization. . . 



B(SX,SY) . CURRENT_COST ;= 0; 
while COUNTER > 0 loop 
COUNTER := 0; 

for I in 6 . . 10 loop 
for J in 1..10 loop 
if I = 6 then 

GRI := B ( I , J ) ; 

DATA_IO. WRITE (OutToMars, GRI) ; 
DATA_IO . READ ( InFromMars , GRI ) ; 



75 



B ( I -1 , J ) := GRI; 

end i f ; 

FIND_MIN ( I , J ) ; 

VOLT := B ( I , J ) . OLD_COST -B ( I , J ) . CURRENT_COST ; 
if VOLT > VOLTA then 

COUNTER := COUNTER +1; 
end if; 

B ( I , J ) . OLD_COST : = B ( I , J ) . CURRENT_COST ; 
end loop; 
end loop; 

integer_inout . put ( counter ) ; text_io . put ( " " ) ; 

GRI . CURRENT_COST := COUNTER; 

DATA_IO. WRITE (OutToMars, GRI); 

DATA_IO . READ ( InFromMars , GRI ) ; 

COUNTER : = GRI . CURRENT_COST; 
integer_inout .put (counter) ; 

text_io . put_line ( " " ) ; 

end loop; 

for I in 1 . . 5 loop 
for J in 1 . . 10 loop 

DATA_IO . READ ( InFromMars , GRI ) ; 

B ( I , J ) := GRI; 

end loop; 
end loop; 



--output of least cost path.. 



loop 

PUSH ( S , B ( GX , GY ) . CURRENT_COST ) ; 
PUSH1 (S1,B(GX, GY) .DIRECTION) ; 
SQ : =B (GX, GY) .DIRECTION; 
if SQ = "NORTH" then 

PUSH ( S , B (GX, GY) .N. DISTANCE) ; 
GX ;= GX+1; 

elsif SQ = "EAST " then 

PUSH ( S , B ( GX , GY ) .E. DISTANCE) ; 
GY := GY+1; 

elsif SQ = "SOUTH" then 

PUSH (S, B (GX,GY) .S. DISTANCE) ; 
GX := GX - 1 ; 

elsif SQ = "WEST " then 

PUSH ( S , B ( GX , GY ) .W. DISTANCE) ; 
GY := GY-1 ; 
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else 
exit ; 
end if; 

exit when GX =SX and GY =SY ; 
end loop; 



DISTANCE 



TEXT_IO . PUT_LINE ( " 

TEXT_IO . PUT_LINE (" 

loop 

POP ( S , E ) ; INTEGER_INOUT . PUT (E) ; TEXT_IO . PUT ( " 
POP ( S , E ) ; INTEGER_INOUT . PUT ( E ) ; TEXT_IO . PUT (" 
P0P1 (SI, El) ; 
if El = "NORTH" then 

TEXT_IO . PUT ( " SOUTH " ) ; TEXT_IO . new_l ine ; 
elsif El = "EAST " then 

TEXT_IO . PUT ( "WEST" ) ; TEXT_IO . new_line ; 
elsif El = "SOUTH" then 

TEXT_IO . PUT ("NORTH") ; TEXT_IO . new_l ine ; 
elsif El = "WEST " then 

TEXT_IO . PUT ( " EAST" ) ; TEXT_IO . new_line ; 
end i f ; 

exit when S. LATEST =0; 
end loop; 



COST DIRECTION"); 



TEXT_IO. CLOSE (INF) ; 
end PRO JO; 
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--Title : PRO Jl. ADA 
--Author : CENGIZ EKIN 
--Date : 20/06/92 
--Revised : 04/10/92 
--Course : THESIS 
--Compiler: MERIDIAN ADA 
--Description : 

with TEXT_IO , COMMON, CHANNELS; 
use COMMON; 
procedure PROJl is 

-- communication channels that are used 

OutToEarth : CHANNELS . CHANNEL_REF := CHANNELS . OUT_PARAMETERS 

( 2 ) ; 

InFromEarth : CHANNELS . CHANNEL_REF := CHANNELS . IN_PARAMETERS 

( 2 ) ; 

B : GRID; GRI : GRID_POINT; 

counter , counterl : integer := 1; 



procedure CAL_WEIGHT (I,J : in INTEGER) is 
begin 

If I = 10 then 

B ( I , J ) .N. WEIGHT := -1; 

B(I, J) .N. DISTANCE := -1; 

B ( I + 1 , J ) . CURRENT_COST :=10000; 
else 

B(I, J) .N. WEIGHT : =1+ ABS ( B ( I , J ) . CURRENT_COST - 
B ( 1+1 , J) . CURRENT_COST ) ; 
end if; 

If J = 10 then 

B ( I , J) .E. WEIGHT := -1; 

B(I,J) .E. DISTANCE := -1; 

B ( I , J+ 1 ) . CURRENT_COST :=10000; 
else 

B(I, J) .E. WEIGHT := 1+ABS (B ( I , J) . CURRENT_COST - 
B ( I , J+l ) . CURRENT_COST ) ; 
end i f ; 

If I = 1 then 

B ( I , J) . S. WEIGHT := -1; 

B ( I, J) .S. DISTANCE := -1; 

B ( I -1 , J) . CURRENT_COST :=10000; 
else 

B(I, J) .S. WEIGHT : =1+ ABS ( B ( I , J ) . CURRENT_COST - B(I- 
1, J) . CURRENT_COST ) ; 



78 



end i f ; 

If J = 1 then 

B ( I , J ) .W. WEIGHT := -1; 

B(I, J) .W. DISTANCE := -1; 

B(I,J-1) . CURRENT_COST :=10000; 
else 

B ( I, J) .W. WEIGHT : =1+ ABS ( B ( I , J ) . CURRENT_COST - B(I,J- 
1) . CURRENT_COST ) ; 
end i f ; 

end CAL_WEIGHT ; 



procedure FIND_MIN (I,J : in INTEGER) is 
begin 

if B(I,J) . CURRENT_COST > ( B ( I + 1 , J ) . CURRENT_COST 

+abs(B(I,J) .N. WEIGHT) ) then 
B ( I , J ) . CURRENT_COST := ( B ( 1+1 , J) . CURRENT_COST 

+abs ( B ( I , J) .N. WEIGHT) ) ; 
B(I,J) .DIRECTION := "NORTH"; 

end i f ; 

if B ( I , J ) . CURRENT_COST > ( B ( I , J+l ) . CURRENT_COST 

+abs ( B ( I , J ) . E . WEIGHT) ) then 
B(I,J) . CURRENT_COST := ( B ( I , J+l ) . CURRENT_COST 

+abs (B (I, J) .E. WEIGHT) ) ; 

B ( I, J) .DIRECTION := "EAST "; 

end if; 

if B ( I , J ) . CURRENT_COST > ( B ( I - 1 , J ) . CURRENT_COST 

+abs (B(I,J) .S. WEIGHT) ) then 
B ( I , J ) . CURRENT_COST := ( B ( I - 1 , J ) . CURRENT_COST 

+abs ( B ( I , J) . S. WEIGHT) ) ; 
B(I,J) .DIRECTION : = "SOUTH"; 

end if; 

if B ( I , J ) . CURRENT_COST > (B(I, J-l) . CURRENT_COST 

+abs (B(I,J) .W. WEIGHT) ) then 
B ( I , J ) . CURRENT_COST := ( B ( I , J- 1 ) . CURRENT_COST 

+abs (B ( I , J) .W. WEIGHT) ) ; 
B(I,J) .DIRECTION := "WEST "; 

end if; 
end FIND_MIN; 



begin 

for I in 1 . . 5 loop 
for J in 1 . . 10 loop 

DATA_IO . READ ( InFromEarth, GRI ) ; 
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B ( I , J ) := GRI; 

end loop; 
end loop; 



--It determines the borders and calculates the weights of the 
edges . . 



for I in 1 . . 5 loop 
for J in 1 . . 10 loop 
C AL_WE IGHT ( I , J ) ; 
end loop; 
end loop; 



--It makes the costs max number in order to use them in 
comparisons for finding 
-- the minimum. . . . 



for I in 1 . . 5 loop 

for J in 1 . . 10 loop 

B ( I , J ) . CURRENT_COST := 10000; 

B ( I , J ) . OLD_COST : = B ( I , J) . CURRENT_COST ; 
end loop ; 
end loop; 



--This part sends dim of matrix, no of 
proccessors , volta , sourceand goal points. 



DATA_IO. READ (InFromEarth, GRI) ; 

N : =GRI . CURRENT_COST ; P : =GRI . OLD_COST ; VOLTA : =GRI . N . WEIGHT ; 

SX : =GRI . E . WEIGHT ; S Y : =GRI . S . WEIGHT ; GX : =GRI . W . WEIGHT ; GY : =GRI . E 
.DISTANCE; 



--cost minimization... 



B (SX, SY) .CURRENT_COST := 0; 
while COUNTER > 0 loop 
COUNTER := 0; 

for I in 1 . . 5 loop 
for J in 1..10 loop 
if I = 5 then 

DATA_IO. READ (InFromEarth, GRI) ; 
B ( 1+1 , J ) := GRI; 

GRI : = B ( I , J ) ; 
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DATA_IO .WRITE (OutToEarth, GRI ) ; 
end i f ; 

FIND_MIN(I, J) ; 

VOLT := B ( I , J ) . OLD_COST -B ( I , J ) . CURRENT_COST ; 
if VOLT > VOLTA then 

COUNTER := COUNTER +1; 
end i f ; 

B(I,J) ,OLD_COST := B ( I , J ) . CURRENT_COST; 
end loop ; 
end loop ; 

DAT A_IO. READ ( InFromEarth, GRI); 
counterl := GRI . CURRENT_COST ; 
if (counter=0) and (counterl=0) then 
GRI . CURRENT_COST 0; 

DATA_IO. WRITE (OutToEarth, GRI); 
else 

COUNTER : =1 ; 

GRI . CURRENT_COST :=1; 

DATA_IO. WRITE (OutToEarth, GRI); 
end if; 
end loop; 



for I in 1 . . 5 loop 
for J in 1 . . 10 loop 
GRI : = B ( I , J ) ; 

DATA_IO. WRITE (OutToEarth, GRI) ; 
end loop; 
end loop; 
end PROJ1 ; 
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# File: makefile 

# "make help" to print option list 

# 

# Complete development cycle: 

# make family -- makes Ada family and library 

directories 

# make -- compiles, links, configures source 

# make run -- run bootable code 

MODE = s 
PROC = 8 

OPTS = / $ (MODE) / 1 $ ( PROC ) 

# make the executable code 

main.btl: mainh . c$ ( PROC )$ (MODE ) proj lh . c$ ( PROC )$ (MODE) 
main .pgm 

@ echo EXPECT 1 WARNING . . . 
iconf /s main. pgm 
@ f:\util\bell 

mainh. c$ ( PROC) $ (MODE) : projO.o pro j Oh . t$ ( PROC )$ (MODE) 
merger . t$ (PROC) $ (MODE) mainh. t$ (PROC) $ (MODE) 
ilink / f main. Ink 

projO.o: common. ada projO.ada 
ada invoke proj 0 . inv, yes 

proj Oh . t$ ( PROC) $ (MODE) : proj0h2.tax projOh.occ 
occam $(OPTS) projOh.occ 

proj0h2.tax: proj0h2.occ 

occam /ta /x proj0h2.occ 

merger . t$ ( PROC) $ (MODE) : merger . occ 
occam $(OPTS) merger. occ 

mainh . t$ ( PROC )$ (MODE ) : mainh. occ 
occam $(OPTS) mainh. occ 

proj lh.c$ (PROC) $ (MODE) : projl.o pro j lh . t$ ( PROC )$ (MODE) 
ilink proj lh . t$ ( PROC) $ (MODE) projl.o adarts8.1ib 
hostio.lib occam8s.lib 
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projl.o: common. ada projl.ada 
ada invoke pro j 1 . inv, yes 



proj lh . t$ ( PROC) $ (MODE) : projlh2.tax projlh.occ 
occam $(OPTS) projlh.occ 



projlh2.tax: projlh2.occ 

occam /ta /x projlh2.occ 

# 

# misc. 

# 

help : 

@ echo Make arguments : 

@ echo make 

@ echo make -n [opt] 

commands 

@ echo make *.o 

@ echo make help 

@ echo make clean 
@ echo make run 

@ echo make check 

@ echo make family- 
directories 

clean : 

del * . ?8? 
del *.tax 
del * . o 
del * . dsc 
del *.btl 

del test_lib\adalib . * 
rd test_lib 
del test_fam\adafam. * 
rd test_fam 

run : 

iserver /sb main.btl 

check : 

check /r 

family : 

ada invoke family . inv, yes 



- make from top level down 

- display but don't execute 

- make Ada object 

- display this list 

- delete all files except source 

- run bootable program 

- check transputer topology 
- make Ada family and library 
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-- File: main.pgm 
# INCLUDE " hostio . inc" 

# INCLUDE " linkaddr . inc " 

PROTOCOL PASS IS INT ; [ 5 ] BYTE : 

#USE "mainh.c8s" 

#USE " pro j lh . c8s " 

CHAN OF PASS Mars2Earth, Earth2Mars : 

CHAN OF SP FromFiler, ToFiler: 

PLACED PAR 

PROCESSOR 0 T8 

PLACE FromFiler AT linkO.in: 

PLACE ToFiler AT linkO.out: 

PLACE Mars2Earth AT link2.in: 

PLACE Earth2Mars AT link2.out: 

[1325000] INT wsl : 

main. harness (FromFiler, ToFiler, Mars2Earth, Earth2Mars, 

wsl ) 

PROCESSOR 1 T8 

PLACE Earth2Mars AT linkO.in: 

PLACE Mars2Earth AT linkO.out: 

[1280000] INT ws2 : 

pro j 1 . harness (Mars2Earth, Earth2Mars, ws2) 
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-- File main . Ink 

-- Purpose: File list for ilink 

mainh . t8s 

merger . t8s 

hostio . lib 

occam8s . lib 

( proj0h.t8s projO.o adarts8 . lib hostio. lib occam8s.lib ) 



-- File: family. inv 

family . new test_fam, overwrite=yes 

lib ( family=test_fam) .new test_lib, overwrite=yes 



-- File: projO.inv 

default . compile library=test_lib 

compile common . ada 

compile pro jO. ada 

def ault . bind library=test_lib, level=bind, warning=no 
bind proj 0 , obj ect = "proj 0 . o " , entry_point= "proj 0 .program" 



-- File: projl.inv 

default . compile library=test_lib 

compile pro jl. ada 

default . bind library=test_lib, level=bind, warning=no 
bind proj 1 , obj ect = "pro j 1 . o" , entry_point = "proj 1 . program" 
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-- File: mainh.occ 
#OPTION "AGNVW" 

#INCLUDE "hostio.inc" 

PROTOCOL PASS IS INT; [ 5 ] BYTE : 

PROC main. harness (CHAN OF SP FromFiler, ToFiler, 

CHAN OF PASS Mars2Earth, Earth2Mars, 
[ ] INT FreeMemory) 



#USE "hostio.lib" 

#USE "proj Oh. t8s " 

#USE "merger. t8s" 

[ 1 ] CHAN OF ANY Debug: 

[2] CHAN OF SP FromAda, To Ada : 

CHAN OF BOOL StopDebug, StopMultiplexor : 
SEQ 

PAR 



-- A multiplexor to combine the debug and normal output, 
so .multiplexor (FromFiler, ToFiler, FromAda, ToAda, 
StopMultiplexor) 

-- A debug channel merger. 

debug. merger ( ToAda [0], FromAda [0], Debug, StopDebug) 

-- A process to invoke the sieve program, 
ws IS FreeMemory: 

SEQ 

proj 0 . harness ( FromAda [1], ToAda [1], Debug [ 0 j , 
Mars2Earth, Earth2Mars, ws ) 

StopDebug ! FALSE 
StopMultiplexor ! FALSE 

so. exit (FromFiler, ToFiler, sps . success) 
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File: merger. occ 



#OPTION "AGNVW" 

# INCLUDE "hostio . inc" 

PROC debug. merger (CHAN OF SP FromFiler, ToFiler, 

[ ] CHAN OF ANY Debug, 

CHAN OF BOOL Stop) 



#USE "hostio. lib" 

-- A debug channel merger and blocker. 

VAL max. debug IS 20: 

VAL number . of . debug IS SIZE Debug: 

INT line. index: 

[256] BYTE line. buffer: 

BYTE value, r: 

BOOL running, reset, s: 

[max . debug] BOOL mask: 

VAL BYTE line. feed IS 10 (BYTE): 

SEQ 

SEQ i = 0 FOR number . of . debug 
mask[i] := TRUE 
running := TRUE 
reset := FALSE 
line. index := 0 
WHILE running 
PRI ALT 

ALT i = 0 FOR number . of . debug 
mask[i] & Debug [i] ? value 

SEQ 
IF 

value = line. feed 
SEQ 

-- Send the complete line, 
so. puts (FromFiler, ToFiler, spid.stdout, 
[line. buffer FROM 0 FOR line. index], r) 
line. index := 0 
mask [i] := FALSE 

reset := TRUE 

TRUE 
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SEQ 

-- Add character to line. 

line . buff er [ line . index] := value 

line. index := line. index + 1 

reset & SKIP 
SEQ 

reset : = FALSE 

SEQ i = 0 FOR number . of . debug 
mask[i] := TRUE 
Stop ? s 

running : = FALSE 
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-- File: projOh.occ 



#OPTION " AGNVW " 

#INCLUDE "hostio.inc" 

PROTOCOL PASS IS INT; [ 5 ] BYTE : 

PROC pro j 0 . harness (CHAN OF SP FromAda, ToAda, 

CHAN OF ANY Debug, 

CHAN OF PASS Mars2Earth, Earth2Mars, 

[ ] INT FreeMemory ) 

# IMPORT " pro j 0h2 . tax" 

[ 1 ] INT dummy . ws : 
wsl IS FreeMemory: 

[3] INT in. program: 

[3] INT out. program: 

SEQ 

-- Set up vector of pointers to channels, 
in . program [ 0 ] := MOSTNEG INT -- not used 

LOAD . INPUT . CHANNEL ( in . program [ 1 ] , ToAda) 

LOAD . INPUT . CHANNEL ( in . program [ 2 ] , Mars2Earth) 

LOAD . OUTPUT . CHANNEL ( out . program [ 0 ] , Debug) 

LOAD . OUTPUT . CHANNEL ( out . program [ 1 ] , FromAda ) 

LOAD. OUTPUT. CHANNEL ( out . program [ 2 ] , Earth2Mars) 

-- Invoke the Ada program. 

-- Assumes the entry point name has been changed to 
"proj 0 .program" . 

proj 0 . program (wsl, in. program, out. program, dummy. ws) 



-- File: PROJOh2.occ 
#OPTION " AEV " 

PROC proj 0 .program ( [ ] INT wsl, in, out, ws2 ) 
[1000] INT d: 

SEQ 

SKIP 
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-- File: projlh.occ 



#OPTION " AGNVW " 

# INCLUDE "hostio.inc" 

PROTOCOL PASS IS INT ; [ 5 ] BYTE : 

PROC proj 1 . harness (CHAN OF PASS Mars2Earth, Earth2Mars, 

[ ] INT FreeMemory) 

# IMPORT "proj lh2 .tax" 

[1]INT dummy. ws : 
wsl IS FreeMemory: 

[3] INT in. program: 

[3] INT out. program: 

SEQ 

-- Set up vector of pointers to channels. 

in . program [ 0 ] := MOSTNEG INT -- not used 

in . program [ 1 ] := MOSTNEG INT -- standard i/o not used 

LOAD . INPUT . CHANNEL ( in . program [ 2 ] , Earth2Mars) 

out . program [ 0 ] := MOSTNEG INT -- standard i/o not used 

out . program [ 1 ] := MOSTNEG INT -- standard i/o not used 

LOAD. OUTPUT. CHANNEL ( out . program [ 2 ] , Mars2Earth) 

-- Invoke the Ada program. 

-- Assumes the entry point name has been changed to 
"proj 1 . program" . 

proj 1 .program (wsl, in. program, out. program, dummy. ws) 



-- File: projlh2.occ 
#OPTION " AEV " 

PROC pro j 1 . program ( [ ] INT wsl, in, out, ws2) 
[100000] INT d: 

SEQ 

SKIP 
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